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ABSTRACT 


Two linearized solvers (time and frequency domain) based on a high resolution 
numerical scheme are presented. The basic approach is to linearize the flux vector by 
expressing it as a sum of a mean and a perturbation. This allows the governing equa- 
tions to be maintained in conservation law form. A key difference between the time 
and frequency domain computations is that the frequency domain computations re- 
quire only one grid block irrespective of the interblade phase angle for which the flow 
is being computed. This is achieved through the use of phase shifted periodic boundary 
conditions. As a result of this and due to the fact that the governing equations for this 
case are steady, frequency domain computations are substantially faster than the cor- 
responding time domain computations. The linearized equations are used to compute 
flows in turbomachinery blade rows (cascades) arising due to blade vibrations. Numer- 
ical solutions are compared to linear theory (where available) and to numerical solu- 
tions of the nonlinear Euler equations. 


CHAPTER I 
INTRODUCTION 


Blade vibration (the aeroelastic problem) and noise generation (the aeroacoustic problem) are 
both undesirable consequences of the unsteady flow processes that occur within a turbomachine and 
are, therefore, important concerns to the designer. The aeroelastic problem has received far more 
attention than the aeroacoustic problem because the former can lead to structural failure of the blad- 
ing and possibly result in extensive damage to the engine [1]. An unsteady aerodynamic analysis in- 
tended for turbomachinery aeroelastic and aeroacoustic applications must be applicable over a wide 
range of structural geometries, mean-operating conditions, and unsteady excitation modes and fre- 
quencies. Because of the large number of parameters involved, the analysis must provide the neces- 
sary unsteady information both efficiently and economically 

To date, three different types of methods have been used to model the unsteady flows 
through cascades [2]. They are : 1) analytical, 2) semi-analytical, and 3) numerical. The degree of 
accuracy of each of these methods depends on how well the mathematical models which represent 
these methods describe the unsteady flowfield for flow through a vibrating cascade. 

The first approach involves solving the partial differential equations which govern the flow in 
a cascade using methods of applied mathematics. Closed form solutions can be obtained for a hand- 
ful of extremely simplified geometries and flow conditions and hence are of limited use for configura- 
tions of practical interest. 

The second approach uses a combination of singularities distributed along the airfoil surfaces 
and wakes to model the steady and unsteady flow through a cascade. This reduces the governing 
equations to a set of integral equations involving kernel functions. Because of the complexity in de- 
termining these kernel functions, the mean flow through the cascade is assumed to be uniform or 
nearly uniform. Hence, the flow model often bears little resemblance to flows found in actual turbo- 
machines. 

The third approach is to discretize the field equations which describe the unsteady flow and 
solve them numerically. The advantage of this approach is that features such as blade loading, arbi- 
trary airfoil shapes, and moving shocks are accounted for. This ability to analyze fairly general ge- 
ometries comes only with considerable computational overhead. The computational methods can be 
further subdivided into 1) methods based on nonlinear governing equations and 2) methods based on 
a set of linearized equations. The methods based on the nonlinear governing equations are almost 
always time-marching, while those based on the linearized equations are time— marching (time-do- 
main) or harmonic (frequency-domain). The nonlinear governing equations are used if the unsteadi- 
ness in the flow cannot be treated as a small-perturbation of the nominal mean flow (steady or peri- 
odic) through the cascade. If the small-perturbation assumption is valid, the nonlinear equations 
can be linearized about the mean flow, resulting in a set of linearized equations for the unsteady per- 
turbations. For most cases of practical interest, the linearized equations are good enough. The har- 
monic method introduces the additional simplifying assumption that the unsteady flow perturbations 
are harmonic in time. This reduces the linearized governing equations to a set of time-independent 
partial differential equations (PDEs) for the complex amplitudes of the unsteady flow properties. The 
time-marching method, on the other hand, makes no such assumptions and solves for the unsteady 
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perturbations by marching the solution in time. Both these approaches have been used in the course 
of the present study. 

This report is divided into seven chapters. Chapter 2 deals with the development of the gov- 
erning equations for the time and frequency domain computations. Development of a high-resolution 
flux formulation and numerical computation of flux Jacobians is presented in Chapter 3. Chapter 4 
covers the solution procedure used while Chapter 5 deals with boundary conditions (BC’s) for the 
time and frequency domain computations. Representative results are shown in Chapter 6. In Chap- 
ter 7, some conclusions are drawn regarding the utility of the methods used in this research. 
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CHAPTER II 

GOVERNING EQUATIONS 


In this chapter, the governing equations for linearized problems (time and frequency domain) 
are developed. The Navier-Stokes equations can most accurately describe the complicated nature of 
steady and unsteady flows through turbomachinery. However, if the Reynolds number is sufficiently 
high and Prandtl number is of order unity, the Euler equations provide a good approximation to the 
behavior of the flow. The nonlinear Euler equations for a dynamic curvilinear grid are obtained by 
transforming the Euler equations from Cartesian coordinates to dynamic curvilinear coordinates 
[3], [4]. These equations form the basis for both the nonlinear (un)steady flows as well as for the linea- 
rization. Section 2.1 develops the linearized Euler equations (time domain) based on the small per- 
turbation approach, the perturbation being about a nonlinear mean flow (obtained from a steady state/ 
unsteady (periodic) solution to the unsteady, nonlinear Euler equations). Section 2.2 develops the 
linearized Euler equations for frequency domain computations. 


2.1 Linearized Euler Equations (time domain) 


For illustration purposes, consider the one-dimensional Euler equations in Cartesian coordi- 
nates, that is, 


dg 

dt 



( 2 . 1 ) 


where q is the Cartesian dependent variable vector and f is the Cartesian flux vector. To start the 
linearization process, the dependent variables are expressed as the sum of two components, namely, a 
meanflow component which, by definition, satisfies the (un)steady nonlinear equations, and an un- 
steady perturbation, i.e., 

?=? + ?' 

where the perturbations q r are assumed small compared to their mean-flow counterparts q. Similar- 
ly, the flux can be viewed as being composed of two components, i.e., 

f = j + r 


where f and f' are as yet undefined, but can be interpreted as mean and perturbed fluxes, respective- 
ly. A definition of f' based on a high resolution flux formulation will be given in Chapter III. 

Using these definitions of q and f, and noting that the mean flow satisfies the equation, 


dq df 

dt dx 


= 0 


Equation (2.1) may be written as 


W , dT 

dt dx 


= 0 


( 2 . 2 ) 


Extending this approach to curvilinear coordinates results in 
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0 


(2.3) 


W | dF 

dr d£ 

£ = £( x , t ) ; Q = J q ; Q = Q + Q' 


F = J(Z x f+Z r q)-, F = F + F' 

This extension is necessary because it is easier to handle arbitrary configurations in curvilinear coor- 
dinates as opposed to Cartesian coordinates. However, F is no longer a pure function of Q; it also 
depends on the metrics, i.e., £ T and £*. Hence, any definition of F f should take into account this 
dependency of F on the metrics. Extending Equation (2.3) to three-dimensions results in 


3Q1 , dF_ + dGl . dK. 

dr * d£ drj dC 


(2.4) 


Equation (2.4) is the linearized Euler equations (time domain) for three-dimensional curvilinear 
coordinates. 


2.2 Linearized Euler Equations (frequency domain) 

Expressing Q f and F' in Equation (2.3) as 
Q' = Re{ Q 0 (£ ) exp(ia>r) } 

F' = Re{F 0 (£)exp(ia>r)} 

results in 

<)F 

0 iu)Q 0 + = 0 (2.5) 

i = 

a) = frequency of vibration 


Q 0 , F 0 = Complex amplitudes of Q' and F* respectively 

Equation (2.5) represents a steady-state equation for the complex amplitude Qq. This equation is 
solved using a time marching method by introducing a pseudo-time derivative , i.e., we solve 


+ {i(o)Q o + if 5 - 0 

Equation (2.6) can be extended to three-dimensions resulting in 


dQ 


+ d(o)Q o + 


dF o 


dG 0 

a| 


+ 


9H o 


= 0 


( 2 . 6 ) 


(2.7) 


Equation (2.7) is the linearized Euler equations (frequency domain) for three-dimensional curvili- 
near coordinates. 
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CHAPTER III 

NUMERICAL FORMULATION 


The governing equations for both the nonlinear and linearized cases can be written in the 

form 


a * ai 


(3.1) 


where d= Q, % = F, (j = G, % = H and y = 0 correspond to the nonlinear case, Q. = Q' , T = F', Q = G' , 
36 = H' and = 0 correspond to the linearized (time domain) case and Q. = Qq, T = Fq, (j = Gq, Dfi = 
Hq and $ = -i(o Qq correspond to the linearized (frequency domain) case. These governing equations 
are hyperbolic in time and can be solved by a time-marching technique. 


3.1 Discretization of the Governing Equations 


A finite volume semi— discretization of Equation (3.1) may be written as 


aa , * , £jj , K * _ 

dt A£ Ay A£ 

where 


(3.2) 


( 3 * (.) = (.),. - (.). , 

2 2 

and i,j, k denotes a finite volume (cell) center. The * ± 1/2 denotes a cell interface and the depen- 

dent variables are assumed to be constant within each computational cell. In the curvilinear trans- 
formation the terms A§, A r] and A£ can be set equal to one and Equation (3.2) can be written as 

^ - -% (3.3) 

at 

where 


% = (<}. 5F + dj q + d k % - S) 


3.2 Flux formulation 

The governing equations are hyperbolic in time, hence any flux formulation process should re- 
spect the direction of information propagation as dictated by the eigenvalues (upwinding). The ap- 
proach used in this research is due to Roe [5] . It is based on a one— dimensional analysis of the Riemann 
problem, the idea being to solve an approximate Riemann problem exactly at each interface. The rest 
of this section is devoted to the development of formulas (first and second order spatial accuracy) for 
the fluxes at the cell faces for both the nonlinear and linearized Euler equations (time & frequency do- 
main). 
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3.1.1 Flux formulas for nonlinear Euler equations 

Details of the development of flux formulas (first and higher order) for the Euler equations are 
available from a variety of sources [5] [6] [7] [8] [9] . In particular, the formulation used in this research 
can be found in [6] . 


3.1.2 Flux formulas for linearized Euler equations (time domain) 


Since the nonlinear unsteady Euler equations were solved using a high-resolution scheme 
based on Roe’s approximate Riemann solver, it is logical to use a similar formulation for the linearized 
equations. The details of the linearization are available in [10]. The final form of the linearized flux 
formulation is: 


5F — Fr Fr — A(q R} q Ly M^) (Qr — Qr) 


(3.4) 


where an overbar Indicates a mean-flow quantity and a prime indicates a perturbed quantity. The 
perturbed flux at any face can then be computed as 


F'u 1 -X + X 


(3.5) 


where 


a.j = strength of the j th wave = • dQ 

t J) = j th eigenvalue of A(q R , Mt) 


t J) , T® = j th left and right eigenvectors of A(q R , q L , MA respectively 
Mt = Vector of metrics = 


Q' = Jq’ + J'q 

dQ = Qr - Ql 

and 2“ denotes summation over the negative wave speeds. As can be seen from Equation (3.5), F' t 

needs to be computed in order to evaluate F. v This can be achieved through a linearization of the 

l + 2 

flux vector, i.e., 

F l = F l + J' A(q L ,M \) q L + J A(q L ,M g ) q' L + J R(q L ,M^)M' ^ 

where 

A(q L > M P = f| I and R(Q l ,M p = ^ | 

Noting that the flux (Fr) can be written as a sum of a mean and a perturbation results in the follow- 
ing definition for 
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F l ' = J' A(q L ,M p q L + J A(q L ,M ^ q' L + J R(q L ,M^)M' ^ 


= A(q L ,M p Q' l + R(Q l ,Ms)M's (3 6 ) 

Similar formulas can be written for computing G' and H\ thus extending this formulation to 
multiple-dimensions. Higher order fluxes are computed in a manner similar to the nonlinear fluxes 
and limiters are used to suppress spurious oscillations. 


3.1.2 Flux formulas for linearized Euler equations (frequency domain) 

Substituting for the “primed” quantities in Equations (3.5) and (3.6), we have 


where 


F ' . i = F 0 exp(icor) and F Q 

l+ 2 


(AQu + RM^ + 


(3.7) 


a- = strength of the j th wave = F • dQ 0 
F = j th eigenvalue of A 


i,r® = j th left and right eigenvectors of A , respectively 
dQ 0 = Q Ro - Q Lo 

Further, 

= M t t exp(icor) 

- r f 

where £ X q etc. are complex numbers and represent the harmonic deformation of the grid due to the 

vibration of the body. Formulas similar to F f can be written for G' and H' also. Higher order accu- 
rate fluxes are computed using simple extrapolation. Limiters are not used because using them 
introduces nonlinearities which prevent certain simplification from being made. 
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CHAPTER IV 

SOLUTION PROCEDURE 


In this chapter, time discretization of Equation (3.2) is carried out. A broad class of difference 
schemes for advancing the solution in time is given by 


+ _ (e _ m , _ 


(4.1) 


Some of the implicit time differencing schemes represented are three-point backward (0=1 ,rp = 1/2), 
backward Euler (0 = l,\p- 0), and trapezoidal (0 = 1/2, xp = 0). Accounting for deforming grids [10] 
results in 


A q n = 


1 At 
1 + ip I J*+ 


T [ (0 - l)& tt - 6'3b n+ 1 ] + ip 


Jn - 1 t /I,?" 71 ”! 

Aa n ~ 1 + a" — 

Jn+1 ^ q Jn + 1 


- q 




(4.2) 


jn+l 

where q is the Cartesian dependent variable vector. Rewriting Equation (4.2) in operator form and 
applying Newton’s method (for time accuracy) to the resulting equation, we have [10] (for the nonlin- 
ear Euler equations) 

'1 +xp I 


where 


0 


-L + 6M + ■ + 6M~ ■ Aq n + l ’ m = - RHS(q n + 1 ’ m ) 
At J 


(4.3) 


RHS(q n + 1 ' m ) = (< ? n+1 ’ m - q n ) - i [(0 - l)R n - 6R n + 1 ’ m ] 


01 


J^Zl Ap n- 1 + 1 ±£ 

0 zlr J n+1 V 0 zlr J ft+1 


1 + ^ 


zlr = 


zlr 

jn + l 


and m represents the Newton iteration parameter Adopting the same approach for the linearized 
Euler equations (time-domain) results in 


"1_ 


+ V J I + dM + • + 6M~ 


, qm + l,m = 


RHS(q' n + 1 ’ m ) 


6 Ax 

where the overbar for the M’s implies mean flow based quantities and 

RHS(q ,n + 1 ’ m ) = q' n+1 > m - q' n ] - J— [ (0 - 1 )R' n - 6R ' n+1 ’ rn ] 


(4.4) 


zlr 1 + xp 


Aq 


tn-1 


8 


The linearized (frequency domain) Euler equations are solved using local time stepping and 
Newton's method need not be used (since we desire only steady-state solutions). The discretized fre- 
quency domain equations can be written as 


1 + ip 

~e~~ 


4~I + dM + 

At 


+ SM 


Aq o = 


(4.5) 


where R 0 is the complex residual. The 6M ± in all these equations can be defined as 

6M + ■ = A + ■ + 6jB + ■ + d k C + ■ 

6M~ ■ = dj A~ ■ + dj B~ ■ + d k C~ • 


4.1 Choice of Flux Jacobians 


The flux Jacobians can be obtained from the flux vector split (FVS) form or the flux difference 
split (FDS) form. The final form of the Jacobians using these two formulations are: 


FVS: 


A* = 


dF ± (q n + 1 ’ m ) 

dgn + l,m 


B' 


dGHq 


n + l,m\ 


dq 


n + \,m 


»± _ 




FDS: 


dF. lX 

= ill 1 

” d qi +P /> 


B : 


dG 

J^2 

dq,- 


r J + k \ 

' J+P s’ 


C ± = 


?aS ) | 

d qk+p/> 


where p = 0 for M + and p = 1 for M and m is the Newton iteration parameter. 

4.2 Solution of Linear Systems 

Equations (4.3), (4.4) and (4.5) can be written in a general form as 

( L+D+U)x=b (4.6) 

where L is strictly lower block triangular, D is block diagonal and U is strictly upper block triangular. 
This system defined by Equation (4.6) can be solved using the “N-Pass” algorithm. The “N-Pass” 
algorithm can be viewed as a relaxation scheme based on the symmetric Gauss-Seidel algorithm: 

( L + D ) x 1 + U x° = b 1 

(D + U)x 2 + Lx 1 =b J lst P ass 

* (4.7) 

( L + D) x 2N_1 + U x w ~ 2 = 61 
( D + U ) x w + L x™-' = b I Nth P ass 
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with x° = 0. The nature of L, D and U can be exploited to reduce the number of matrix-vector multi- 
plications to be carried out [10]. 
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CHAPTER V 

BOUNDARY CONDITIONS 


The problem of determining a solution to a partial differential equation when both initial and 
boundary data are present is called an initial-boundary value problem [11]. The solution to the Euler 
equations on finite domains falls under this category of problems. Boundary conditions are vital in de- 
termining the success or failure of a CFD code. Incorrect boundary conditions can often lead to wrong 
solutions or at least drastically alter convergence rates. Boundary conditions (BC’s) can be classified 
into three major categories: 1) Far field BC’s (inflow/outflow) 2) Solid wall BC’s (no flow across a surface) 
3) Block-to-block BC’s (transfer boundary conditions). 

The Euler equations are hyperbolic in time. This means that the characteristics determine the 
direction of propagation of information from or to the boundaries. Therefore, it is only natural to use 
a formulation based on characteristic variables to determine the physical and numerical boundary 
conditions to be imposed at any given boundary. The formulation and implementation of far field condi- 
tions for steady flows is identical to that used in [12] and will not be repeated here. 

This chapter is divided into four sections. The first section deals with the development of non- 
reflecting far field boundary conditions for unsteady flows (time domain) as well as for the frequency 
domain. Sections 2 and 3 deal with development of impermeable wall BC’s for nonlinear and linearized 
flows respectively. The last section explains the multi-block strategy for the unsteady (time domain) 
cases and how the computations can be performed in just one grid block for the frequency domain cases. 

5.1 Non-Reflecting Boundary Conditions 


The two-dimensional Euler equations (Cartesian coordinates) in nonconservative, quasilinear 
form can be written as 


dU , dU , z. dU _ n 

-rr + a r b ~~z = U 

dt dx dy 


(5.1) 


where U = ( q , u, v,p) T and 



U Q 0 O' 


— 'i 

c 

o 

'O 

o 


0 U 0 i 


0 V 0 0 

a — 

Q 

0 0 u 0 

b = 

00 » i 


0 yp 0 u 


0 0 yp v 
< > 


= 0 


Equation (5.1) can be linearized to yield (for small perturbations): 

dir 

dt 

where 

U' = \jq',u',v',p'] t 


, n dU' h dU' 

-r a —7 r 0 — r — 

dx dy 


(5.2) 


and a and b are constant coefficient matrices based on the steady flow variables. Non-dimensional- 
ization of Equation (5.2) results in: 



'u 1 0 O' 


'v 0 1 O' 

a = 

0 u 0 1 
0 0 u 0 

b = 

0 v 0 0 
0 0 v 1 


010b 

k J 


001b 


A Fourier analysis of Equation (5.2) (as described in [13]) yields a set of (four) eigenvectors repre- 
senting an entropy wave, a vorticity wave, an upstream running pressure wave and a downstream 
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running pressure wave. Based on this analysis, the characteristic variables for unsteady flows can be 
written in terms of the perturbation variables as: 


' C 1 


1 0 or 


V 

c 2 


0 010 


du 

c 3 


0 101 


dv 

c 4 

* 


0 - 1 0 1 


3? 


where dg, du> dv and dp are the perturbations from the mean flow and the c; ’s are the characteristic 
variables. Using these eigenvectors and assuming locally one-dimensional flow at the boundary, the 
boundary condition, expressed in terms of the characteristic variables is simply 

C/ i=0 


for all n corresponding to incoming waves. This means that, for subsonic inflow, the amplitudes of 
incoming unsteady characteristics (cj, C 2 , C 3 : assuming co-directional inflow) are set to zero and the 
outgoing characteristic is computed using Equation (5.3). For subsonic outflow, c 4 = 0 (again, assum- 
ing co-directional outflow) and the remaining characteristics are computed. Once the characteristics 
are known at theinflow as well as outflow boundaries, the perturbations can be determined using an 
inverse transform. 

The primitive variables can now be calculated because the perturbations from the mean flow 
(i.e., dg , du , dv and dp) are known. Once the primitive variables have been determined, the conserved 
variables can be calculated. The strategy for the frequency domain case is very similar to that described 
here because the governing equations can be cast in terms of the nonconserved quantities, thereby 
yielding a set of equations identical to (5.2) except that the primed quantities are replaced by the corre- 
sponding complex quantities. 


5.2 Solid Wall Conditions for Nonlinear Euler equations 


Casting the non-conservative form of the nonlinear, Euler equations into characteristic vari- 
able form and assuming locally one-dimensional flow results in 


dr 


+ A 1 


Wk 
’■ dk 


= 0 


k = £ ,rj or £ 


(5.4) 


where are the characteristic variables. The details of the derivation and implementation are 
available from a variety of sources including [4], [12], [14] etc. 

5.3 Solid Wall Conditions for Linearized Euler Equations 


The linearized Euler equations can be written as 


dQ' 

dr 


+ I 

k =|,?7 andf 


= 0 


(5.5) 


where an overbar indicates a mean quantity and a prime a perturbation (a real number for time do- 
main and a complex number for frequency domain) and K = A for k = §, K = B for k = r\ and K-C for 
k = £, A, B and C being the flux Jacobians in the £, rj and £ directions respectively Casting Equation 
(5.5) into characteristic variable form and assuming locally one-dimensional flow yields 


where 


dWk_ 

dr 


+ At 


Wjk 
’■ dk 


= 0 


(5.6) 
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and 


W k = T k,0 Q' = l w 'k,l> w 'k,2> w 'k,3> w 'k,4>w' k ,5l 


K = T k A k T ;> 

The subscript 0 (zero) indicates a reference condition. An impermeable surface is characterized 
by zero normal velocity. Therefore, the first three eigenvalues are zero, the fourth positive and the fifth 
negative. Using this in Equation (5.6) yields 


(w' k>1 )b = ( w' kl )r 

(5.7) 

( w ' k$)b ~ ( W 'k,2 )r 

(5.8) 

( w 'k,3\ = ^ W 'k,3^r 

(5.9) 

]k x U + k y V + k z W + k t j b = 0 

(5.10) 

( w 'k;*)b = ( w 'k,d r 

(5.11) 

or 


( w ‘ k, b ~ ^ w 'k,5^ r 

(5.12) 


where the subscripts b and r denote boundary and reference values respectively. In this case, the 
reference value is chosen as the center of the first cell off the boundary. If the point r is in the positive 
k direction from the boundary, Equation (5.12) is used, otherwise Equation (5.11) is used. Equations 
(5.7) through (5.9) are imposed conditions, i.e., they help in obtaining a set of five simultaneous equa- 
tions that can be solved. Note that even though Equation (5.10) is written in terms of total quanti- 
ties, a perturbed form of the same is used in the solution process. This is because the rest of the 
equations are in perturbed form. This system of equations (real or complex perturbations, but with 
real coefficients) is solved for the perturbations at the boundary. 

5.4 Block-to-Block Boundary Conditions 

Figure 5.1 represents a single grid block used for calculations with steady flows and flows with 
in-phase blade motions. All the cases considered in this study involved cascades. The boundaries that 
extend upstream and downstream of the airfoils are called periodic boundaries . Therefore, a periodic 
boundary condition is imposed between lines A-B and E-F, and C-D and G-H. 



Figure 5.1 Periodic boundaries in a typical cascade 

For harmonic blade motions with constant phase difference (interblade phase angle, a) between 
adjacent blades, multiple blade rows ( N of them) are required to simulate the flow field. The number 
of blade passages required for arbitrary a is given by N = smallest integer { ( 360/a ), 360/(360 - a) }. 
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For example, for a = 180°, two grid blocks are required. For the frequency domain case, multiple grid 
blocks are unnecessary because a phase shifted boundary condition can be applied at the periodic 
boundaries thereby enabling computations to be performed in just one grid block. 
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CHAPTER VI 
RESULTS 


Various test cases were used to validate the applicability of the present method for solving 
problems involving aeroelastic and aeroacoustic phenomena. The test cases selected cover a wide 
range of geometries as well as flow conditions. This chapter is divided into three main sections. The 
first section deals with uniform mean flows, while the second one deals with nonuniform mean flows 
(flows with blade loading). The last section discusses the performance of the linearized solver with 
respect to the nonlinear solver. In both cases, the unsteadiness in the flow field is due to blade mo- 
tion. The blades oscillate according to: 

h y — hty sin (cot + o) for tangential (plunging) motion 

a = a 0 sin(a>£ + a) for torsional (pitching) motion 

where ho y and cto are the respective amplitudes, co the circular frequency of oscillation and a the inter- 
blade phase angle- These oscillations are captured using a locally deforming grid technique [10] . 

As a part of the code validation process, linearized results are compared with linear theory for 
a cascade of flat plates. Linearized Euler computations are also compared with nonlinear Euler solu- 
tions for a cascade of airfoils, for which no theoretical results are available. The computations typically 
required 4 cycles of motion with 200 time steps/cycle (for time domain) and 200 pseudo time steps (for 
frequency domain) to reach a periodic/steady state. A schematic of a typical two-dimensional cascade 
defining various terms used in this chapter is shown in Figure 7.1. 




a = Angle of attack 


Figure 7.1 Schematic of a typical cascade 


7.1 Uniform Mean Flows 


For the set of cases considered in this section, the mean flow is uniform, i.e., the cascade under 
consideration is made up of infinitely thin flat plates and is at zero degree angle of attack. For these 
cases, theoretical results ([15], [16]) for supersonic as well as subsonic axial flows are available, hence, 
they form a basis for the validation of the present method. The computer program of [15] is based on 
the theory developed by Lane [17] . Unless otherwise mentioned, all computations were performed on 
a 71x51 H-grid, with 50 points on either face of the flat plate. 


15 


7.1.1 High solidity cascade, Plunging motion, M 00 =2.61 


Linearized results (time domain) are first compared to available theoretical results for a high- 
solidity cascade with a supersonic axial flow. The cascade is oscillating in plunge, has a gap-to-chord 
ratio of 0.311 and a stagger of 28°. The freestream Mach number is 2.61, the interblade phase angle 
is 90°, and the reduced frequency (k) based on semi-chord is 0.5. The maximum CFL used was 5. As 
can be seen from Figure 7.2, excellent agreement is seen between theory and the linearized solver (time 
domain). The location as well as the strength of the real and imaginary parts of the discontinuities are 
captured accurately by the linearized solver. 



Distance Along Chord x/c 

Figure 7.2 Unsteady pressure difference on a flat plate cascade due to plunging 
motion (M 00 =2.61 / s/c=0.311, 7=28°, k=0.5 and a = 90?) 

7.1.2 Low solidity cascade, Pitching motion, M 00 =0.8 

This case involves a cascade oscillating in pitch at a reduced frequency of 1.5 with an interblade 
phase angle of 270°. The cascade geometry is identical to that used in the previous test case and the 
free stream Mach number is 0.8. A higher resolution grid (121x41) was used in this computation and 
the maximum CFL was 3.25. The comparison is shown in Figure 7.3. The computed results (time do- 
main) are in good agreement with the theoretical predictions. 

7.2 Nonuniform Mean Flows 

Nonuniformity in a mean flow field can be caused either by the cascade being at an angle of at- 
tack or by the cascade being composed of airfoils (not infinitely thin flat plates) or a combination of the 
two. Since the mean flow is not uniform any more, it needs to be computed. This is done using the high 
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Figure 7.3 Unsteady pressure difference on a flat plate cascade due to pitching 
motion (M^O.8, s/c=1.0, y=45°, k=1.5, o = 270° and pitching axis 

@ 30% chord) 

resolution, flux difference split, time-marching solver, which upon convergence yields the required 
steady state/unsteady (periodic) mean flow field. All steady flow computations were performed on a 
single block grid, while the unsteady calculations required multiple grid blocks depending on the inter- 
blade phase angle. The frequency domain computations were performed om a single block grid. This 
rest of this section deals with the so called 10th standard configuration. 

The 10th standard configuration (10th std. cfg.) consists of a cascade of modified NACA 0006 
airfoils. The gap-to-chord ratio is 1.0 and the stagger angle is 45°. All computations were performed 
using a 121x41 H-grid, with 60 points each on the pressure and suction surfaces. The maximum CFL 
used in the steady mean flow computations was 5000 (with local time stepping), while for the unsteady 
cases it was 85. The cascade geometry and the multi-block grid used in the computation is shown in 
Figure 7.4. 


7.2.1 Tenth standard configuration: Subsonic case 


The cascade is in a flow with M oo =0.7 and at a 10° angle of attack. The mean flow through the 
cascade is subsonic and reaches a maximum Mach number of 0.92 on the suction surface. A back pres- 
sure (ratio of p ex it to ptotal) of 0.87 was used to achieve the required Mach number (0.7 in this case) at 
the inlet. It was also observed that the back pressure required was a function of the configuration being 
considered as well as the grid used for the computation. The steady pressure distribution for this case 
is shown in Figure 7.5. For the unsteady calculations, the cascade is set into motion in plunge or pitch. 
The amplitudes of plunge and pitch are 1 % chord and 1° respectively. These cases are considered next. 

Plunging motion, k = 0.6435. o = -90°: 

The cascade is oscillated in plunge at a semi-chord reduced frequency of 0.6435 and an inter- 
blade phase angle of -90°. The real and imaginary parts of the first harmonic of the unsteady pressure 
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distribution are shown in Figure 7.6 and Figure 7.7 respectively. As can be seen from these two figures, 
the agreement between the linearized and nonlinear solution is very good. 



Figure 7.6 Real part of unsteady pressure distribution on subsonic 10th std. 

cfg. due to plunging motion (k=0.6435, o=-90°) 


7.2.2 Tenth standard configuration: Transonic case 

The second case is a transonic cascade (M^O.S, 0=13°). This case is characterized by a fairly 
strong in-passage shock. A back pressure of 0.87 was found to be sufficient to achieve the required 
Mach number of 0.8 at the inlet. The pressure distribution for the transonic case is shown in Figure 7.8. 
The unsteady cases are considered next. 


Plunging motion, k = 0.6435. a = -90°: 

The cascade is pitching at a reduced frequency of 0.6435 and an interblade phase angle of -90°. 
Figure 7.9 shows the real part of the unsteady pressure distribution, while Figure 7.10 shows the imag- 
inary part of the same. As can be seen from these two figures, reasonably good agreement is achieved 
between the nonlinear and linearized Euler solutions. It was observed that for the cases considered, 
the linearized solutions were scalable with respect to the amplitudes of pitch or plunge as the case may 
be, but the nonlinear solutions exhibited a dependency on the amplitudes used (mainly near disconti- 
nuities). 


6.4 Comparison of the Various Approaches Used 

The linearized (time-domain) code required about 14 minutes of CPU time (for a 4 block, 121x41 
grid) on a CRAY-YMP to reach a periodic state. For the same problem, the nonlinear code required about 
20 minutes of CPU time. The linearized (frequency domain) code, on the other hand, required only 
about 1 minute of CPU time (Note that computation in the frequency domain case requires only 1 grid 
block). 
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Imaginary Unsteady Pressure 
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Figure 7.9 Real part of unsteady pressure distribution on transonic 10th std. 

cfg. due to plunging motion (k=0.6435, o=-90°) 
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Figure 7.10 Imaginary part of unsteady pressure distribution on transonic 
10th std. cfg. due to plunging motion (k=0.6435, o=-90°) 




CHAPTER VII 
CONCLUSIONS 


Two linearized Euler solvers (time and frequency domain) based on a high-resolution flux 
difference split scheme have been presented. The linearized equations are maintained in conserva- 
tion law form and the same solution procedure is used to solve the nonlinear as well as the linearized 
equations. Newton subiterations helped improve time-accuracy of both the nonlinear as well as the 
linearized (time domain) solutions. 

The results presented in Chapter 6 prove the feasibility of using the present approaches in solv- 
ing small amplitude aeroelastic and aeroacoustic problems. Though it has not been tested, no practical 
difficulties are anticipated in extending this method to three-dimensional problems. The present meth- 
od performs well over all flow regimes, i.e., subsonic, transonic and supersonic, and various cascade 
configurations. 

As anticipated, the frequency domain based method yields significant performance improve- 
ment without any degradation in the quality of the solution, thus making it a valuable design tool. 
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